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ABSTRACT 

The linear stability of compressible plane Couette flow is investigated. The correct and 
proper basic velocity and temperature distributions are perturbed by a small amplitude nor- 
mal mode disturbance. The full small amplitude disturbance equations are solved numeri- 
cally at finite Reynolds numbers, and the inviscid limit of these equations is then investigated 
in some detail. It is found that instability can occur, although the stability characteristics of 
the flow are quite different from unbounded flows. The effects of viscosity are also calculated, 
asymptotically, and shown to have a stabilizing role in all the cases investigated. Exceptional 
regimes to the problem occur when the wavespeed of the disturbances approaches the ve- 
locity of either of the walls, and these regimes are also analyzed in some detail. Finally, the 
effect of imposing radiation-type boundary conditions on the upper (moving) wall (in place 
of impermeability) is investigated, and shown to yield results common to both bounded and 
unbounded flows. 


iThis research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS1-18605 while the authors were in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1. Introduction 


Incompressible plane Couette flow possesses, perhaps, the simplest exact solution of the 
Navier-Stokes equations, and, (probably as a consequence) the study of the stability of this 
flow has been the subject of considerable attention over the years. Numerical studies of the 
linear stability problem have been carried out by Grohne (1954), Gallagher and Mercer (1962, 
1964), Deardorff (1963), Davey (1973), and Gallagher (1974); however, all these studies found 
no evidence of instability. A number of analytic studies have also been carried out on this 
problem. Wasow (1953) showed that the flow is stable at all streamwise numbers (a) if the 
Reynolds number (Re) is sufficiently large. The stability of the flow at low Reynolds number 
was demonstrated by Synge (1938). Dikii (1964) proved that all modes with wavespeeds equal 
to the average of the wall velocities were stable and indeed that the imaginary component 
of the wavespeed Im {c} < — a/Re. The first general proof of stability appears to be due to 
Romanov (1973). He showed that all normal modes of the linear problem are damped for 
a > 0, Re > 0. Exact solutions of the Orr-Sommerfeld equations have been obtained by 
Reid (1979). 

Another issue that has been studied is the question of a continuous spectrum. Case 
(1960, 1961) showed that the time dependent inviscid problem has a continuous spectrum 
which decays in time as 1/t, this spectrum arising as a direct consequence of the singularity 
of the inviscid equations at the critical layer. Shivamoggi (1982) presented an example of a 
continuous spectrum which decayed as an exponential in time. 

A detailed analysis of modes for the large Reynolds number limit for general mean flows is 
given by Morawetz (1952), which has implications for plane Couette flow. She classified the 
modes into three sets. In the first set are eigenvalues which approach the inviscid eigenvalues 
in the limit of infinite Reynolds number. Since there are no discrete inviscid eigenvalues in 
the case of plane Couette flow, this first set is empty. In the second set, there is an infinity 
of eigensolutions unrelated to the inviscid problem, satisfying 

| c — c n | < A(aRe) - 2, 

where A is a constant, and c n is a root of 

ry 2 i rnr 

,m L [•tvs')-*-] 

where n is any integer, U 0 (y) is the velocity profile, a is the wavenumber, and the flow 
extends from y — y\ to y — y 2 * These modes are always stable (if Re is sufficiently large). 
In the third class are eigensolutions of the viscous problem that, as orRe * oo, approach a 
finite V-shaped strip in c-space defined by one branch of 

Rl / W1 [i(Uo{y) — c))*dy = 0 (1.3) 

Jyc 

and one branch of 

Rl I [i{U 0 (y) - c)]*dy = 0, (1-4) 

where U 0 (y c ) = Rl{c}. These modes are either neutral or stable. 


( 1 . 1 ) 

( 1 . 2 ) 
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However, there does exist a dichotomy between the theoretical and computational results 
described above and experimental results at large Reynolds numbers, in which instabilities 
are certainly observed (Taylor 1936, Reichardt 1956, Robertson 1959). A number of attempts 
have been made to explain this through the inclusion of nonlinear terms. Investigations of 
this type include the work of Watson (1960), Eckhaus (1965), Hains (1967), Kuwahara 
(1967), Reynolds and Potter (1967), Ellingsen, Gjevik, and Palm (1970), and Lessen and 
Cheifetz (1975). These studies used a variety of techniques, including varying degrees of 
mathematical rigor, and led to a number of conclusions (some partly contradictory), although 
on balance, the evidence was that finite amplitude effects could, indeed, destabilize the flow. 
Additionally, Romanov (1973) considered the nonlinear initial- value problem and showed 
that there is a unique solution which is asymptotically stable if the norm of the initial 
disturbances in Sobolev space is sufficiently small. 

The linear stability of compressible flows is considerably less well understood than corre- 
sponding incompressible flows. Most of the work (which has for the most part been based on 
the parallel mean flow assumption) has been with regard to boundary-layer flows (e.g., Lees 
and Lin 1946, Mack 1963, 1965a, b, 1969, 1984, 1987, 1990) and more recently to jets and 
shear layers (e.g., Tam and Hu 19S9a,b, Greenough et al. 1988, Papageorgiou 1990, Jackson 
and Grosch 1989, Mack 1990). 

The first (of a number) of distinguishing features of the stability of compressible flows 
was found by Lees and Lin (1946). This relates to the replacement of the (classical) inflec- 
tion point condition in the streamwise velocity profile (for the existence of neutral inviscid 
modes), by the generalized inflection point (G.I.P.) condition which involves the mean den- 
sity distribution also. This condition relates to the existence of neutral “inviscid” modes (i.e., 
modes of wavelength comparable to the characteristic scales of the mean flow, for example 
the displacement thickness in the case of boundary-layer flows, or lateral extent in the case 
of confined flows). 

Another effect of compressibility is the frequent existence of multiple modes of instability 
in supersonic flows. In addition to the modes related to the G.I.P., a further class of modes 
may exist; these are linked to the so-called supersonic neutral modes (as opposed to subsonic 
neutral modes associated with G.I.P.’s). Such neutral modes are readily distinguishable, 
because of the oscillatory nature of the eigensolutions in the free stream, while subsonic 
neutral modes exhibit exponential decay of their eigensolutions in the free stream. Mack 
(1963, 1965a, b, 1969, 1984, 1987, 1990) has shown how such modes are extremely important 
in compressible boundary-layer flows, particularly in the case of boundary layers on cooled 
walls, where these modes are often associated with the largest amplification rates. Such 
modes have also been shown to exist and be important in other classes of flow by Tam 
and Hu (1989), Greenough (1989), Zhuang, Kubota, and Dimotakis (1990) (confined two- 
dimensional supersonic mixing layers), Macaraeg (1990), and by Jackson and Grosch (1989) 
(compressible mixing layers). However, to date, a simple mechanism that explains why and 
precisely under what conditions these types of mode exist has yet to be propounded. Some 
general remarks concerning stability at large Reynolds numbers have been made, however, 
by Morawetz (1959), in particular, the presence of inner viscous regions and their nature 
was vigorously established in both the incompressible and compressible cases. In the case 
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of stable disturbances, with finite damping, it was shown that a finite viscous region in 
the interior of the fluid, no matter how large the Reynolds number, exists. In the case 
of unstable disturbances, with finite amplification, there is no inner viscous region if the 
Reynolds number of the mean flow is large enough. In the case of neutral disturbances (for 
constant Prandtl number) it was also rigorously shown that if a viscous region does exist, 
its width shrinks to zero with an increase in Reynolds number. 

The object of this study is to analyze some of the characteristics of compressible plane 
Couette flow. Although analytic expressions for the mean streamwise velocity f/ 0 (y ) and 
temperature T 0 (y) are not available for general viscosity laws, T 0 (y) may be expressed as 
a second-order polynomial in Uo(y), and hence it is relatively straightforward to generate 
mean flow profiles, under different conditions. Some work has previously been carried out 
investigating the compressible stability of plane Couette flow. Glatzel (1989) has considered 
this problem, and has carefully studied the least stable eigenmodes. He also derived critical 
Reynolds numbers for the instabilities that were found. However, this work is for the special 
case of constant viscosities, density and pressure (implying constant temperature also), which 
lead to a fourth-order system (rather than a sixth-order system). Girard (1988) considered 
the same problem, and also assumed constant viscosities (although he did allow for variations 
in density and temperature); instabilities were again found to occur. In this paper, we choose 
to consider the stability of the mean flow profile described above, which is a correct and 
proper solution of the full compressible Navier-Stokes equations. We show that the details 
of the mean flow profile have a profound effect on the stability of the flow. 

The layout of the paper is as follows. In Section 2 we formulate the problem, stating 
our fundamental assumptions and equations of motion. In Section 3 we derive the equations 
of motion for the basic flow, in which the problem is reduced to that of a straightforward, 
although nonlinear, first order system that may be solved by means of standard numerical 
means. In Section 4 we derive the (full) small amplitude disturbance equations, neglecting 
only terms in perturbation amplitude squared, and we describe a numerical scheme to treat 
this system, together with a number of numerical results. In Section 5 we consider the 
inviscid limit of these equations; we show how the so-called “generalized inflection point” is 
relevant in this context, and under what circumstances we can expect such a point to occur 
with our basic profile. 

Numerical results for the inviscid problem are presented in Section 6, and, most impor- 
tantly, it is shown that unstable modes are possible, although according to inviscid theory, 
there are many regions where many modes are neutrally stable. These neutral modes are 
investigated further in Section 7, in which the effects of viscosity are included, and shown 
to always have a stabilizing role for high Reynolds numbers. However, this study also raises 
important questions regarding the applicability of our results, in particular in regions of 
changeover from the neutral to non-neutral state. This region is investigated in some detail 
in Section 8. Because of the apparent discrepancy in results between bounded and un- 
bounded stability analyses, in Section 9 we consider a change in boundary conditions on the 
upper (moving) wall, from one of impermeability to one of radiation. Finally, in Section 10 
we present some conclusions. 
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2. Formulation 


We assume that we have a compressible, Newtonian, perfect fluid between two infinite 
parallel planes defined by y* = 0 and y* = h. The x* axis is taken to lie in the plane of the 
lower wall. We take the fluid to have density p*, viscosity p*, second coefficient of viscosity 
£* (which may be taken to be zero for a monotonic gas). The upper wall has velocity Uqo 
(a subscript oo refers to unperturbed conditions on the upper wall), taken parallel to the 
plane of the wall, while the lower wall is at rest. Although the basic flow will be taken to 
be one-dimensional, later we shall consider two-dimensional perturbations of this flow. The 
velocity components are taken to be u* = (it*, u*) in the x * and y * directions respectively, and 
the pressure and temperature are written as p* and T* respectively. We non-dimensionalize 
velocities with respect to f/oo , lengthscales with respect to h, density with respect to p^, 
viscosities with respect to p^, temperature with respect to and pressure with respect 
to where the gas constant R* = C v - C„, and C p and C v are the specific heats 

at constant pressure and volume respectively; non-dimensional quantities are denoted using 
the same notation as corresponding dimensional quantities, except without the superscript 
asterisk. 

The continuity equation may then be written 


| + £ ( H + ^(H = °. (2-0 

Here, and throughout the paper, we assume that the appropriate dimensional timescale is 

0(h/U oo). 

The momentum equations are written 


du du du 

m +a ai + v B^ 


1 dp 

7 


+ 


1 f d 


Re 1 dx 


du 

2 *Tx + AV ■ U 


dv dv dv 

Tt +u Tx + v Ty 


a_ 

+ dy 
1 dp 


( du dtAl 1 

** + dx) J ’ 


iMldy 


1 / d 


^ dv 

r> l a — 1" ^ ' u 

Re [dy [ dy 

Here the Reynolds number Re is defined by 


d [ ( du dv \ 1 1 

+ di V^ + dx)\ J ' 


Re = UooP^h^ 
foo 


the Mach number M ^ by 


Mo o = VoJyJlR'T^ 


and the ratio of specific heats by 7. 


( 2 . 2 ) 


(2.3) 


(2.4) 

(2.5) 
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The non-dimensional form of the energy equation used is 


'dT dT dT\ ( 7 - 

b u b v — - - — 

dt dx dy I V 7 


dp dp dp 

di + u di + v d^ 



'pdT 
a dy_ 


, d_ 

dx 


pdT 
a dx 


My - 1)M* 
Re 




1 / du chA 2 A / du du \ 2 

+ 2 \dy + dx ) + 2p \<9x ^ dy) 


Here cr is the Prandtl number 

a = p*C p /K*, 


( 2 . 6 ) 

(2.7) 


where K* is the coefficient of heat conductivity. The equation of state is simply 


P = PT. (2.8) 

We assume that viscosity depends solely on temperature, and in particular we assume Suther- 
land’s law 

^ = r3/2 (tTc) ’ (2 ' 9) 

where C is a constant. Finally, A = ( — 2/3 fi and the Stokes assumption ( = 0 is assumed 
throughout the paper. 

In the following section we consider the basic flow which we expect to depend on y only. 


3, Compressible Couette Flow 

We seek a solution of (2.1) - (2.6) which is dependent on y only, together with a constant 
mean pressure. By continuity we must have v = 0. We then seek a solution of the form 


u 


= U 0 (y), T = T 0 (y), p = p 0 (y). 


(3.1) 


We then have 


HoToy 


{^oUoy)y 0 , 

+ (7 - = 0, 


(3.2) 

(3.3) 


(3.4) 


subject to 

£/o(0)=0, f/ 0 (l) = 1, 

T o {0) = T w , T 0 (l) = l. 

It follows immediately from (3.2) that the shear stress r is a constant through the profile, 

i.e., 

r = poUay = constant. (3-5) 

The energy equation may then be written as 


p-Toy + (y-l)MlU 0 r 

L cr 


0. 


(3.6) 
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Defining the recovery temperature (i.e., the wall temperature with adiabatic conditions) 


and the recovery factor 



(3.6) may be integrated and written in the form 


To = T r (r + (1 — r)U Q - (l - ^-) . 


(3.7) 

(3.8) 

(3.9) 


Equations (3.5) and (3.9), together with an appropriate viscosity law (we used (2.9)) and 
boundary conditions (3.4) then completely determine the problem; r is unknown a priori, 
and must therefore be determined as part of an iterative process. The mean pressure profile 
is constant. 


4. The Small-Amplitude Disturbance Equations 


Here we take the solution to be that of Couette flow (as discussed in the previous section) 
together with a small amplitude perturbation. More specifically, we write 

U = U 0 (y) + 6u(y)E + O(6 2 ), 

V = 8av(y)E + 0(8 2 ), 
p = 1 + 8p(y)E + 0(8 2 ), 

P = Po(y) + 6p(y)E + 0(6 2 ), (4.1) 

P = Po(y) + 8p(y)E + 0(8 2 ), 

C = Co (y) + 6((y)E + 0(8 2 ), 

T = T 0 (y) + 8T(y)E + 0(8 2 ), 

where 

E = exp [ia(x — ct)} , (4-2) 

and 8 is the (small) amplitude of the perturbation. If we take terms 0(8 ) in (2.1) - (2.3), 
(2.6) we obtain the following full, small disturbance equations: 


— icp + ipou + iUop + vpQ y + poVy — 0, 


po [— iacu + iaUou + otvUo y ] + 



p = 


(4.3) 


— — ^ — 2poa 2 U + Aq —Ct 2 U + ioiVy + Poy [u y + iotv\ + Po + lQ; fiy] + PyUoy + pUoyy 


po f— ia 2 v + m 2 f7ofi] + 


(4.4) 


Py 


7 ML 


^2/Xoy^y T 2/iofiyy 4“ -^Oy (lO.lL T ^y) -^0 {lOLlLy T Vyy ) + Po ( iau y — a 2 if) + iapUo y }, 

(A 
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(4.6) 


p 0 \-iacf + iaUoT + avT 0y \ - [~ioicp + iaU 0 p] = 

{^ 0 Tyy + POyTy + P^yy + PyT[)y ~ 0 7 } 

+ — — ^PoUoyUy + p(Uoy) 2 + 2poiotUoyV^ . 

The perturbation equation of state is then 

P = Top + (4.7) 

No-slip and impermeability conditions are applied to the velocity perturbations at both walls, 
namely 

«(0) = «(1) = 6(0) = 6(1) = 0. (4.8) 

The temperature boundary conditions are T(l) = 0, while on the lower wall, 

T (0) = 0 (4.9) 

for a heated/cooled surface and 

r)T 

—— = 0 on y = 0 (4-10) 

ay 

for an insulated wall. 

We now present results based on the full set of viscous compressible linearized equations, 
i.e. (4.3)-(4.7) above. 

The stability results are obtained from a spectral temporal linear stability code written 
by Herbert (1990). The code assumes a global representation of all variables in appropriate 
basis functions which can vary from variable to variable. The basis functions for the velocity 
components satisfy the imposed boundary conditions. Velocity perturbations are set to zero 
at y = 0 and y = 1. Thus, the velocity basis functions are linear combinations of Chebyshev 
polynomials. In terms of the Chebyshev polynomials T n (y), these basis functions are defined 
as 

U n (y) = T n+ i(y)-T n (y). (4.11) 

Boundary conditions for the density and temperature perturbations waves are of more general 
nature, so these variables are represented as a standard Chebyshev series. In this section, 
only adiabatic results are shown. Thus, the y derivative of the temperature perturbation is 
zero at the lower wall. The upper wall is insulated; therefore the temperature perturbation 
is set to zero. The density at the walls is obtained by integration of the continuity equation. 

The numerical algorithm departs from standard collocation schemes which operate in 
physical space. Taking the perturbation density (p) as an example, let 

n 


7 


and its first derivative 


(4.13) 


d f>{y) _ dT n (y) 

dy dy 

be evaluated at each collocation point. Inserting p and the other variables into the linearized 
stability equations leads to a system of linear equations for the coefficients of the basis func- 
tions. From these, one can either calculate a global spectrum using routines from the IMSL 
library or obtain a single eigenvalue at a sequence of parameter values using a continuation 
technique. Details are given in Herbert (1990). 

Figures 1-5 all correspond to = 2. The dependence of the eigenvalues of inviscid 
character is shown in Figures 1 a-c which correspond to a = 0.05,0.10 and 0.20 respectively, 
all at Re = 2 x 10 5 , and obtained with 100 collocation points. It is clear from the sequence of 
figures that one inviscid eigenvalue is approximately stationary, while the remaining inviscid 
eigenvalues originate from very large values of Rl{c} at very low a. A magnification of figure 
lc (shown in figure Id) shows the region of the spectrum corresponding to 0 < < 1. 

One notes the general “Y”shape of the spectrum as discussed by Morawetz (1952). The 
spectrum is composed of three “Y” shaped pieces. These are due to the structure of the 
viscous terms of the two momentum equations and the temperature equation. The coupling 
between the equations and the non unit Prandtl number leads to the non superposition of the 
“Y” singular curves as Re — » oo. Before drawing any general conclusions however, we show 
spectra at Re = 7 x 10 s in Figure 2, leaving all other parameters unchanged with respect to 
Figure 1. As expected, higher resolution is required to properly resolve the spectrum. The 
sequence of Figures 2a-d correspond to resolutions of 75, 100, 125 and 150 collocation points 
respectively. Although Morawetz (1952) predicts that there is a set of “viscous” eigenvalues 
which lie on the edges of the “V” part of the “Y”, packed with a density proportional to 
\A*Re, it is clear that at M = 2, the triple point of the U Y” is cut off by two horizontal 
bands of eigenvalues. The word horizontal is used here as a qualitative description. At 
the lower resolution oi N = 75, the bands are slightly wider then for higher N . It is not 
completely clear whether these bands will remain or disappear as TV — ► oo. A definitive 
answer would require further investigation into the properties of the viscous component of 
the spectrum at high Reynods numbers and high Mach numbers. Comparing Figures 2b and 
2d, it is clear that one effect of insufficient resolution is the splitting of the spectrum from 
the vertical part of the Y (which is indicated by “|”). However, such a splitting is not 
always indicative of loss of resolution (see Figure 5 at a = 0.1, M x = 5, Re = 2 x 10 6 ). At 
a = 3.5, Re = 7 x 10 5 , the horizontal bands become more evident (Figure 3). Note that the 
density of eigenvalues along this band has increased significantly. Figure 4 shows a blow up 
of the upper left hand corner of the stable diagram. The scales along Rl{c } and Im{c } are 
identical to bring out the real shape of one of the “triple points”. At TV = 150, one notices 
the left branch of the “Y” which is at the theoretical angle of 60°. As TV changes from 125 
to 150 nodes, the horizontal bands shift slightly in the direction of increased stability. Thus, 
the final converged state of these bands is inconclusive for the chosen set of parameters. 

We now turn our attention to a similar set of spectra at M 0 0 = 5 which seem to exhibit 
a more complex character than their M <*, = 2 counterparts. At Re = 2 x 10 6 , a = 0.1, the 
three “Y” shaped curves are clearly present. These structures subsist in an unmodified form 
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(i.e., without the horizontal bands) at higher Reynolds numbers at the higher M w . One 
“Y” is approximately symmetrical about the Rl{c} = 0.5 axis. The “|” structure has split 
into three pieces, one along Rl{c) = 0.5, and the other two curves placed symmetrically 
about Rl{c } = 0.5. The actual locus of points along these curves shows a similarity with 
the continuous temporal boundary-layer spectra discussed in Ashpis and Erlebacher (1989). 
As the Reynolds number is increased to 5 x 10 6 , the resolution requirements simultaneously 
increase. Only when the vertical locus of eigenvalues along Rl{c } = 0.5 appears, is the 
resolution adequate, at least for the part of the spectra above it, although this statement is 
not quite true near the horizontal bands. The unchanging vertical position of these bands 
as N increases from 125 to 150 indicates that its presence is not an artifact of a loss of 
resolution. 

To complete the picture, Figures 7 and 8 show the spectra at M <*, = 5, a = 3.5 at 
Re = 2 x 10 5 and Re = 7 x 10 5 . In both cases some inviscid modes (those near the 7m{c} 
axis) have moved into the 0 < Rl{c } < 1 range. As each mode crosses Rl{c } = 0 or 
Rl{c } = 1 regions, a critical layer develops; this is examined in Section 6. Once again, as 
the Reynolds number increases, the width of the horizontal bands increases, and resolution 
studies indicate that they do not disappear as N — ► oo. 

In anticipation of the results presented in the following sections, we plot the evolution 
of the mode I phase velocity (finite Rl{c} as a — ► 0) as a function of a. This is carried 
out for several resolutions and at two different Reynolds numbers. Although inviscid theory 
suggests stronger instabilities at M = 5, viscous calculations indicate that the resolution 
requirements become much more severe. We therefore restrict ourselves to M 0 0 = 2. Both 
at Re = 7 x 10 5 and Re = 1.4 x 10 6 . Im{c } exhibits a local maximum near a = 4. This 
mode becomes less stable at higher Re, while a shifts to a slightly lower value. As expected, 
a higher resolution is required at Re = 1.4 x 10 6 (N = 150), as opposed to N — 100 at 
Re = 7 x 10 5 . Unfortunately, we are not able to increase Re beyond 1.4 x 10 6 and still obtain 
accurate results at resolutions not exceeding N = 150. Figures 9 and 10 show the variation 
of Rlci and Rlc r for both Reynolds numbers. The phase velocity has been scaled by \/Re to 
bring out the one-half power scaling law. Further discussion of Figure 9 is deferred until the 
asymptotic theory at large Re has been presented. 

In the following several sections, we examine solutions of the stability problem in the 
limit Re — > oo. Guided by a number of previous viscous theories, together with the work of 
Morawetz (1952, 1959), we expect, in general, the solution to develop in the form 

u = u 0 (y) + Re _ = ui(y) + 0(Re _1 ), 

v = v 0 {y) + Re" 2 {ii(j/) + 0(Re -1 ), 

f = To(y) + Re _ 2j’i(y) + 0(Re _1 ), (4-14) 

p = Po{y ) + Re * P\(y) + 0(Re ), 
c = Co + Re" 2 Ci + 0(Re -1 ). 
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5. Inviscid Disturbance Equations 


Taking equations (4.3) - (4.7), and then retaining lowest order terms in Re (see Mack 
1984, for example), we obtain the following two first-order equations for the leading order 
normal velocity (u = u 0 ) and pressure ( p = p 0 ) perturbations respectively 


VQy 


UpyVQ _ ipo {To — — Cp) 2 } 

Uq - Co 7 Uo- c o 


POy 


ia 


( Uq - Cp)u 0 , 


iMl T 0 

where the appropriate boundary conditions to be applied to this system are 

vo(0) = 5 0 (1) = 0, 


(5.1) 

(5.2) 


(5.3) 

(5.4) 


implying 

Po v (0) =Po y (l) = 0, 

where Co is the complex wave-speed of this system. Equations (5.1) and (5.2) may be 
combined to yield the single second order equation (see for example Mack 1984, 1987) 


d j (bp — co)i>oy — UoyVp ^ ^ {Uq — cp) 


dy \ To - Ml(U 0 - coY j 


To 


-Up. 


(5.5) 


for 5p. 

Before investigating any of the above systems numerically, it is interesting to study 
the significance of so-called “generalized inflection points” which are found to be highly 
significant in the context of compressible flows. If we multiply (5.5) by v£ (the complex 
conjugate of u 0 ) and then subtract from the resulting equation its complex conjugate, the 
following equation is obtained: 


d i 

\ {Uo — Co)5o y UoyVo) 

i _ d J 

\ {Uo — Cq)Vq v — UoyVy ] 

U 0 - Cody) 

l X J 

[ U 0 - Cq dy | 

l x* j 


where 


X — To — Mlo(Uo ~ co) 2 , 


(5.6) 


(5.7) 


and a superscript asterisk here denotes a complex conjugate. After some algebra, this may 
be written 


(5.8) 

(5.9) 


S -A 

VOy 

„ d 

O 

Voy 

-vv'l 1 d 

\Uoy] 

1 d 

\U 0 y 

°iy 

. X . 

dy 

.x*. 

° i To - Co dy 

. X . 

ZX, 

*<? 

1 

1 

. x* 


Writing 

c 0 = Rl{co) -f i Im{co] 

then the neutral state corresponds to 7m{cp} — ► 0. In this limit, (5.8) may be written 
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' VpVOy ~ ^Oy 
X 

Using arguments similar to those of Duck (1990), then (i) as /m{co} — ► 0, the left-hand-side 
of the equation is zero except possibly at y, (where Uo(yi) = RI{cq}), (ii) the term in the 
parentheses on the left-hand-side of (5.10) must be zero at both y = 0 and y = 1, while (iii) 
the right-hand-side is clearly non-zero unless 


y i 


voviVi) 2 th 
Uoy(yi) 


C¥) 


(5.10) 


v> 


d_ 

dy 


u. 


Oy 


LTo \ v= 


= 0. 


y=Vi 


(5.11) 


Thus, in order to avoid an inconsistency, we must have (5.11). 

Most importantly, there does exist a difference between the present (bounded) flow con- 
figuration, and that of unbounded flows, in that in the present situation, (5.11) is necessary 
if 

0 < i?/{co} < 1, (5.12) 

i.e., a critical layer must exist inside the flow, while in the case of unbounded flows condition 
(5.11) only holds if the wavespeed is “subsonic”, i.e., 1 — 1/M ^ < 7?/{c 0 } <C 1 4" 1/Moo. This 
is because (5.11) is a direct consequence of the zero velocity perturbations at the domain 
boundaries. There is absolutely no distinction made here between supersonic and subsonic 
modes. Thus any neutral inviscid mode satisfying (5.12) must be associated with a general- 
ized inflection point. The above says nothing about neutral inviscid modes outside the range 
of (5.12). It is also worth noting that authors who implement Dirichlet boundary conditions 
(in place of radiation boundary conditions) on the disturbance terms in truncated infinite 
domains, may well experience difficulties computing non-inflectional supersonic modes, since 
the arguments above suggest that a generalized inflection point is necessary for supersonic 
disturbances, a condition which is clearly erroneous in the unbounded case. 

Utilizing (3.5) and (3.9) in (5.11) yields, 


f r °» T 
i/OT 


[»otTo + ft 0 ] 


= 0 . 

y=y% 


(5,13) 


If the term inside the square brackets is zero, then by Sutherland’s law (2.9), we must 
have 

[3T 0 + 5 C) y=yt =0, (5.14) 

which is clearly inadmissible. Consequently, the only way that a generalized inflection point 
will occur is if the mean temperature profile has a local extremum. If we invoke Sutherland’s 
law and (3.9), then we either require 


U 0 y(y .) = 0, 


which is clearly not possible on account of (3.5), or 


Uo(Vi) = 


(1 ~ rJTr 

2 [T r ~ 1] ' 


(5.15) 


(5.16) 
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Since T r > 1, and the flow is unidirectional, then this condition cannot be satisfied unless 


0 < U 0 (yi) < 1, 


which implies 


1 - Ml 
1 + 


< r < 1, 


(implying that the lower wall must be cooled below adiabatic conditions). 


(5.17) 

(5.18) 


6. Inviscid Disturbance Results 

For this set of computations, the basic profile for velocity and temperature, namely Uo(y) 
and T 0 (y) was determined by means of a fourth-order Runge-Kutta scheme. The eigenvalue 
problem posed in (5.1) - (5.4) was solved similarly, with Newton iteration being used to 
update the complex wavespeed Cq so that all boundary conditions were satisfied. 

In a number of computations (specifically those for which 0 < RI{cq } < 1 ), it was found 
necessary to deform the integration contour into the complex y plane (in particular below 
the real y axis), in a manner described by Lees and Lin (1946). This was undertaken by two 
independent ways. The first, as suggested by Mack (1965b), involves obtaining the mean flow 
profile (U 0 (y),To(y)) along the real y axis, and then using Taylor series expansions to obtain 
the mean flow profile off the real y axis; this detour is made close to where Uo{y ) = Co- The 
second approach (which was generally used in preference) obtains the mean flow solution 
itself in the complex y plane, thereby eliminating the errors associated with truncation of 
the Taylor series. Specifically, the mean flow was obtained for 0 < y < y x (i.e., on the real y 
axis), yi < y < y\ — iy 2 , y\ — iyi < y < y$ — iy 2 , yz — iy 2 < y < J/a and thereafter back along 
the real axis jfa < y < 1 . J/i,?/ 2 ,J /3 were all taken to be real and positive (although this is 
not essential) and were chosen to avoid the computation proceeding too close to the critical 
layer. Comparison of results using the two approaches proved a useful check of the accuracy 
of our results. A further check on our results was that in addition to solving (5.1), (5.2) we 
also solved the adjoint system (see (7.29), (7.30) below), and also (5.5). 

The first results we present is for the case M ^ = 2 and adiabatic lower wall boundary 
conditions (here and in all cases we took 7 = 1.4, a — 0.72 together with the Sutherland 
constant C = 0.5). Results for Rl{co} are shown in Figure 11 . The results display a number 
of interesting features. There appear to be many modes, which, however, can be divided 
into two distinct families, the first corresponding to 72/{co} > 1 as a — * 0. All these modes, 
with the exception of one (mode I) have Rl{co} — > 00 as a — ► 0; we refer to these as the 
upper family. These modes all suffer a monotonic decrease in RI{cq} as a increases, and 
the results suggest that some finite value is approached in the limit a —*■ 00 . The second 
family is defined by RI{cq) < 0 as a — » 0, and the f?/{co} all increase monotonically with an 
increase in a, and eventually become positive. We refer to these modes as the lower family 
of solutions. Again, all except the first of these modes (mode II) appear to be unbounded 
as a — » 0 . 

A further, important feature is that in all cases for which RI{cq} > 1 or Rl{c 0 } < 0, all 
these modes are neutrally stable to within machine accuracy, i.e., 7m{co} = 0. However, 
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since there can be no generalized inflection point for the mean flow under consideration, then 
there can be no neutral modes with 0 < Co < 1. Indeed, it is found that in the case of the 
upper family of solutions, once Rl{co} < 1, then these modes become stable according to 
our inviscid calculations, while in the case of the lower family of solutions, once R1{cq } > 0, 
then these modes become unstable according to our inviscid calculations; the location at 
which Rl{c o} either drops below 1 or rises above 0 is marked on Figure 11 by a circle. The 
distribution of 7m{co} for mode I is shown in Figure 12 (other upper family modes also 
have negative values of 7m{co}, but many orders of magnitude smaller). The distribution of 
7m{co} for mode II is shown in Figure 13 (other modes of the lower family have considerably 
smaller values of 7m{c 0 }, although these are also positive). Next we go on to consider the 
aforementioned trends as a — > 0 suggested by our numerical results. 

The existence and behavior of modes I and II as a — ► 0 is easy to confirm. If we set 
a = 0 in (5.5), and integrate once, we obtain 

(Co - co)v 0y - UoyVo = I<[To - MUUo - Co) 2 ], (6.1) 


where K is an arbitrary constant. This equation may be integrated once more to yield 


Vo 


K (Uo-Co)J* 


[T 0 - MjjUo - 
(Co - co) 2 



( 6 . 2 ) 


If this is to satisfy the boundary condition on y = 1 (the boundary condition on y = 0 is 
already satisfied by (6.2)), we must then have 


l 


0 (Co -Co) 


-,dy = K 


(6.3) 


This relationship was solved for Co (for given Co(i/), T 0 (y), and M 0 c ) and confirmed the 
numerical results as a — > 0. Indeed, if we take a model problem, i.e., 


UoW = V, T 0 (y) = 1 (6.4) 

then the solutions of (6.3) are 

1 ± [1 + 4 /Ml)'/* 

co= . (6.5) 


These modes will be discussed further in Section 7. 

The behavior of modes III and higher as a — ► 0 is different, but nonetheless straightfor- 
ward to confirm. In these cases we have Co — ► oo as a — ► 0, and a balancing of terms in (5.5) 
in this limit demands 


Co — ot I coo + O(l), (6*6) 

Vo = uoo + 0(a). (6*7) 


The equation for vqq is then 


VQyy + 


U 00 


Mlv 


^00 


To 


= 0. 


(6.8) 
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If, to make further progress analytically we assume |coo| > 1, implying we are focusing 
attention on the higher nodes (alternatively we could assume the model problem, with T 0 = 
1), then the WKB solution is 

voo = ATJ | e icooM °° So r o _1/2 ^ - e -^oMoo /; T~ 1 / 2 d y | ? ( 6 9 ) 

where A is a constant. Here the boundary condition on y = 0 has been applied, while the 
other boundary condition on y = 1 demands 


Ax) — 


mr 

Voo/cX^V 


( 6 . 10 ) 


where n is any (large) positive or negative integer. This expression clearly illustrates the 
multiplicity of modes. 

The second set of results we show is for the higher Mach number case = 5, with 
the adiabatic boundary conditions on the lower wall retained. Figure 14 shows a number of 
results for Rl{c 0 }. Although these seem much the same as the corresponding distributions 
for Mqo = 2 , there are some differences, in particular with mode I a, as defined in Figure 
14, which initially corresponds to mode I in Figure 12 . However, unlike the corresponding 
results for M^ = 2 , it turns out that there exist other (stable) modes in the vicinity of mode 
I A- Some of these modes are shown in Figure 15 (on a magnified scale) and are denoted by 
Ib and Ic ■ These modes were difficult to generate numerically in regions of a where Rl{c 0 } 
was very close to unity and where the decay rates were quite small; however, the existence 
of these modes was confirmed using the three different formulations of the inviscid stability 
problem described above, together with the condition ( 6 . 3 ). Indeed, it is quite likely that 
other such modes exist, although the present investigation failed to yield any. However, it 
turns out that these modes are all stable, as shown in Figure 16. Note that modes I a and 
Ic have, in places, values of Co exceedingly close to each other, but nonetheless distinct; this 
compounded the difficulty associated with these modes. 

Modes II, and higher, on the other hand did exhibit the same qualitative behavior as the 
corresponding M ^ = 2 results. In particular, mode II becomes unstable at a ~ 1.85, with 
the distribution of 7m {co} shown in Figure 17. Other higher modes of this lower family are 
also unstable, but with substantially smaller growth rates. 

The final set of results in this section relates to the case M^ = 5, as previously, but with 
the (cooled- wall) boundary condition 7o(0) = 1. It is worth noting that even in this case, the 
uniform temperature together with a linear velocity profile is still not a proper solution of 
the governing equations, except in the very special limit as 7 — » 1 . Nonetheless, the solution 
does have some important properties, namely that 


To(y) = 7i(l - y), 


( 6 . 11 ) 


together with 

U Q {y) = 1 - U 0 (l - y). ( 6 . 12 ) 

As a consequence of these symmetries, the mean flow has a single G.I.P. at y = (with 
Uq = |). These properties turn out to have interesting implications for the inviscid 
stability of the profile. 
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Figure 18 shows the variations of RI{cq } for the first eight modes. As in the previous cases, 
there again exist two distinct families, with all except two of the modes having the property 
_> oo as a -> 0. Furthermore, all modes are neutrally stable if RI{cq } < 0 or 
Rl{c 0 } > 1. Significantly, due to the aforementioned symmetries about y = the behavior 
of RI{cq] for the upper family of modes is merely the mirror image of the corresponding 
lower family member about y = These symmetries also yield the result that /m{co} is 
precisely the same for corresponding modes at the same wavenumber. Figure 19 shows the 
distribution of /?n{co} for modes I and II. These modes begin as neutrally stable, then as 
i?/{c 0 } drops below unity/rises above zero, both modes become stable. As a increases still 
further, both modes have /?/{co} — * | which is reached at a w 5.5. Im{co } then becomes 
positive, implying unstable modes. This is entirely consistent with our G.I.P. arguments, 
which predict neutral modes having Co = Thereafter, 7?/{co} for mode I increases, and 
decreases for mode II. 

Figure 20 shows the distribution of Im{co} for modes III and IV, which over the range 
shown are seen to be stable. It is possible that regions of instability exist for higher a, 
although the computations become increasingly difficult as a — > oo. 

In the following section, we go on to consider regions where viscous effects are likely to 
become important, particularly in determining the stability properties of the flow. 

7. The Viscous Correction 

The results described in the previous sections show that there is a family of solutions of 
(5.1) - (5.4) comprising neutrally stable modes over a wide range of a, which have either 
Co > 1, or Co < 0, and hence have no critical layer. The question then arises as to the effect of 
viscosity on these modes - whether it plays a stabilizing or destabilizing role in this problem. 
We investigate this question next. 

The viscous correction to the problem arises due to the thin layers that occur as a result 
of the violation of the no-slip condition (and also the temperature condition) on the walls 
y = 0 and y = 1. These conditions are 


u = v — 0 

on y = 0 and y — 1, 

(7.1) 

T = 0 on y = 1, while on y = 0, 

^ = 0 
oy 

(7.2) 

if the wall is insulated, and 

II 

o 

(7.3) 

if it is heated or cooled. 

Now from (4.4) and (4.6) we see that 

on y = 0 we have 


u 0 

T o (0)po(0) 

IMS,* ’ 

(7.4) 

and x 

M 

To(0)p(0), 

(7.5) 
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while on y = 1 we have 


= ToWoO) 

0 tM^co-1)’ 


(7.6) 


and 

ft = f 1 ^-) Io(l)#b(l). (7.7) 

Thus, generally these expressions fail to satisfy the appropriate wall conditions described 
above. 

In order to remedy this, we require a (Stokes-like) layer of thickness 0(Re - *) on both 
walls. Considering first the layer on y = 0, defining 


Y = yRe* = 0(1), 

(7.8) 

together with 


u = U(Y) + 0( Re"*), 


v = Re-*R(F) + 0(Re- 1 ), 


f = f(Y) + 0( Re-*), 

(7.9) 

P = P(Y) + 0( Re-*), 


p = }(Y) + 0( Re-i), 


then taking 0(1) terms in (4.3) - (4.7), we have 


—iacoU ~ iaP 

ft(0) = I ' lUy¥ ~ IMI’ 

(7.10) 

Py = o, 

(7.11) 

!Co ' ,+ T„( 0) + T(0) H ' =0> 

(7.12) 

T Tyy+( ) tacoPo -0, 

T o (0) (7 V 7 / 

(7.13) 

P — To{0)p + j> o(0) > 

(7.14) 

where jil — /^o(0). On account of (7.11) the pressure must be constant across 
so we write 

this layer, and 

P = po(0) = p L . 

(7.15) 

The solution to (7.10) is then 


jj - h c (*^51^) Y \ 

I [ 

(7.16) 
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In the case of an insulated wall on y = 0, the temperature perturbation is given by (7.5) 
across the layer (i.e., it is constant), while for the heated or cooled wall conditions, using 
(7.13) we obtain 

^~ l ^p L T{ 0){l-e“(' r ^) V 


T = 


(7.17) 


7 


Using the continuity equation, (7.12), we obtain 


V = icoPLY - icopL 


' 7-1 ' 

. 7 


r + 




ipLl o(0) 

' iM^co 


<Y + 


-icxcp 

7o(0 


[ —iacpa 1 2 

MlT 0 (O)J 

.-(wSr) 1 _ i 


(7.18) 


where b = 0 in the case of an insulated wall, while b - 1 for a cooled/heated wall. The 
corresponding layer on the upper wall is quite similar. Defining the lengthscale 

Y = (y- l)Re», 

together with (7.9), then the results for Y = 0(1) may be simply inferred from those for 
Y - 0(1) by (i) replacing l Y' by ‘-Y\ (ii) replacing L pY by ‘pi/’ = p(l), (iii) replacing ‘co’ 
by “co - 1” and (iv) replacing p L and T o (0) by unity. Finally the (key) result for V is then 
given by 

V = i(co-l)puY 


-i(co - 1 )pu 


' l ~ 1 

v 7 . 


Y - 


-ia(co - 1)<t] 2 


.[-iorfco-lppy _ 1 


ipu 


7 M^(c 0-1) 


Y 


We shall later require the behavior of (7.18) and (7.19) as Y 
respectively. We have 

7 — 1\ 


V 


Y ^ icopL - icqPl 


(7.19) 


1 LHc.(co- 1)1*K ill 

-ia(co - 1)]* J j 


oo and Y 


-co, 


Pl 


7^Co. 


icoPLb (^) 

I ^ 

— iacpo 2 

t 0 (0) ml . 


ipL'I o(0) 


7M^co 


(7.20) 
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and 


V 


— ► v 0L Y + v 1L as y -+ oo, 
l v ; [-ic<(c 0-1W* 


■M 


7M^(co-l)[-iO!(co 


“oil} 


»(co-l)pt/( V) ^ 


CoJ7^ + V XU as Y 


-oo. 


We now seek the leading-order viscous correction term to (5.1), (5.2) for y = 
using the expansions (4.14). The variables v\ and p\ turn out to be determined 
inhomogeneous form of (5.1) and (5.2), namely 


V\y - 


ip i 


UpyVl 

Uo — Co Uo — Co 

ia 2 7 MI 


To - Ml{Uo - cof 


where 


Ply + 
Ri = 


7 Ml 

(Uo - Co)Cl = CiJ?2, 

, ipoTo ip 0 

(U 0 - Co) 2 7 Ml(U 0 - co) 2 7 ’ 

Mia 2 Vo 




T 0 

UoyVo 


Ro = 


To 


The boundary conditions on this system are also inhomogeneous, namely 

hi(0) == v XL , 

f»i(l) = v w , 


(7.21) 

(7.22) 

(7.23) 

0(1), 

by an 

(7.24) 

(7.25) 

(7.26) 

(7.27) 

(7.28) 


(the boundary conditions for p ly may be obtained from (7.25), although this is not necessary 
for the following), where V\i and V\u are defined by (17.21) and (7.23) respectively. 

To obtain ci, we use the condition of solvability of the system (7.24) - (7.28). For this 
we require the adjoint to the system; if we denote v + and p + as the adjoint functions, then 
these are to be determined by 


+ , CV+ , ia 2 (Uo - cJ-fMl + 
vJ + — 1 " 


Uo — Co 


To 


0, 


Pt 


n T 


Uo — Co 


T^-MUJJo-cof 

7 Ml 


= 0. 


The boundary conditions to be applied to this system are that 


p + (0) =p + (l) = 0. 


(7.29) 

(7.30) 


(7.31) 


Ci is then given by 


^(l)^ - U + (0)f?iL 
ti[R l v+ + R 2 p+]dy' 


(7.32) 
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where v+ and are the complex conjugates to and respectively. 

Given the nature of the solution for 72, and R 2 when Cq is real, it is easy 

to show that we must have 

Ci = Im{ci}(±l + i), (7.33) 

where the positive sign is taken for Cq > 0, the negative sign is taken for cq < 0; note 
that 7m{ci} > 0 for instability. However, in all the computations we performed, without 
exception, we found 7 m {ci} < 0. Results for Im{ci] for the adiabatic case with == 2 
(modes I and II) are shown in Figure 21. The distribution becomes unbounded as a — » 0, 
and also as a — > a 0 , the point at which Co — > 1 or 0; these distributions are typical of 
all modes. Figure 22 shows the corresponding distributions for the adiabatic case with 
Moo = 5, and exhibit, qualitatively, the same behavior. (The corresponding computation for 
Moo = 5, T o (0) = 1 was also carried out, and yielded qualitatively the same picture.) We 
shall defer discussion of the limit a — ► a 0 until the following section, but let us now consider 
the limit as a — ► 0, of modes I and II from which it appears that our expansions cease to be 
uniformly valid. This is most clearly seen by the apparently singular behavior of the viscous 
correction to the complex wavespeed as a -> 0, together with the 0(a~ 2 ) growth of the 
wall-layer thickness in this limit (see (7.10) for example). Thus when a = (^(Re^ 1 ) the wall 
layers will fill the entire channel, and the disturbances become purely viscous in nature. 


Specifically, if we write 


a = Re 

(7.34) 

then the governing equations (4.3) - (4.6) reduce to 


-icp + ip 0 u + p 0 v + iU 0 p + vpoy = 0, 

(7.35) 

zap 

ap 0 [ icu + iU 0 u + vU 0 y] + 

1 1 CO 

(7.36) 

— f^O^yy ~f"~ PO'j'TU^y ~f~ ft'Q'J'^Qy^y y 


o' 

II 

(7.37) 

ap 0 [—icT + iU 0 T + vT 0y ] ^ ^ ^ [iotU 0 iac]p 


POrP , POT m rh , TpoTrr , TyPor 

= lyy + ±0yly + 1 0yy J t 

(7.38) 


a a a <7 

-\-T Pot tTq v + T pori') — l)-VL^oy> 

where it has been implicitly assumed that viscosity is a function of temperature only. The 
problem then reduces to a basically viscous system. 

However, for modes III and higher, Co = O(o -1 ) as a — + 0, and since the Stokes layer 
thickness remains 0(R~ 2 ), it does not fill the entire channel. Additionally, as a — ► 0, 
7m{ci} = 0(cv _1 ), and so the expansions (4.14) remain valid. 

Returning to Figure 9, we note the general excellent agreement between the values of 
\/Re Im{c) obtained numerically, and /m{ci} obtained from the above asymptotic theory, 
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for the case = 2 (mode II). The values of c r obtained using the two approaches are 
almost indistinguishable. We do note a deviation between the results as Co — ► 0, caused 
by the (expected) breakdown in the asymptotic theory in this limit (and as cq — > 1). This 
aspect is investigated in the following section. 

8. The Nature of the Solution as Co — ♦ 0 (or Co — ♦ 1) 

We consider here the nature of the complex wavespeed as cq — > 0 (the results for Co — + 1 
may be simply inferred from those of Co — ♦ 0). Our previous results indicate that as Co — ► 0, 
(with a — > a 0 ) then (i) c x , the viscous correction to the wavespeed becomes unbounded, 
(ii) the leading order wavespeed cq becomes non-neutral as a increases above a 0 , and (iii) 
the thickness of the wall layer increases. Thus this region is likely to be a regime of some 
interest, to which we now focus our attention. 

More specifically, guided by our numerical results, if we write 


a = o-o + ec*i, (8.1) 

where 

£ = Re - s, (8.2) 

then we expect 

C = £Ci+0(£2). (8.3) 

The solution in the “core,” i.e., away from the wall layers is then expected to develop in the 
form 

v = 6o(y) + eMy) + ' ’ ‘ , (8-4) 

P = Po(y) + epi{y) + ■ • ■ • (8.5) 

The leading order system is then 


v 0y - 


U , OyVO 


ip o 

1MIU 0 




(8,6) 


or symbolically 
and 

or symbolically 
together with 


£i{u 0 ,Po} = 0, 


Poy . • 2 Uovo n 

WE + ,a °W = °’ 


Aftio.Po) = 0, 


«o(0) = v„(l) = 0. 


(8.7) 

( 8 . 8 ) 

(8.9) 

( 8 . 10 ) 


Notice that (8.6) and (8.8) imply that po = 0(y 3 ) as y — t 0; this system effectively determines 
the value(s) of a 0 for which Cq = 0. Turning to the next order system, we find 


U 0 Ci 


= Uoy-TJT + 


*P 0 


US iMIUS 


^1 [To + M^Uq] — c 1 R 1 , 


(8.11) 
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C2{vi,Pl} = 


ialcivo 2ia 0 aiU 0 Vo 


T 0 To 

The boundary conditions to be applied to this system are 


— C\ /?2 + OL\ R3. 


*i(l) = 0, v 1 (0) = Ap 1 (0), 


( 8 . 12 ) 


(8.13) 


where A is to be determined later. Using conditions of solvability, we must have that 


» _ -t) + (0)Api(0) - Qq /q 1 RoP^dy 

“ ~ " /„'[£,«+ + R,p*\iy 


(8.14) 


where p + (y), v+(y) are the adjoint functions to (8.6), (8.8), namely those determined by 
(7.29), (7.30) with a — a 0 , cq = 0. 

We now consider the effect of the wall layers. The upper wall layer remains of thickness 
0(Re“2) and as such plays no role to this order. The interesting changes are related to the 
lower wall layer, where now we must have 

Y = e~ l y = 0{ 1) (8.15) 

as the crucial scale, wherein to leading order 

(u,v,p) = (Uo,eV 0 ,eP 0 ). (8.16) 


Taking the leading order terms in the continuity and momentum equations we find 

iUo T Voy = (8.17) 

Par = 0, (8.18) 

i°to[Uo{0 )Y — c\]Uq + q 0 Cq(0)Vo = Uqyy ~ 1 Q o / 0 > (8.19) 

with boundary conditions 

r) o (0) = v o (0) = 0, 

and as Y — * 00. 

Uo -* -iv ^(0), (8.20) 

where the latter condition arises from a proper matching with the y = 0(1) solution. If we 
now differentiate (8.19) with respect to Y , invoke (8.17) and (8.18) we find 

U Q yyy ~ ia 0 [U^( 0)Y ~ C,]U Q y = 0, (8.21) 


whose solution is 

Urt = C l Ai{[ia 0 UttO)]HY -<*)}, 

where C\ is a constant to be determined, and C\ = Ci/f/o(0). 

Utilizing the boundary condition on Y = 0, and also evaluating (8.19) on Y = 0 leads to 


taohjl ^{[*ao^(0)]?(y-cQ}dr 
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Vo may now be determined by continuity (8.12) together with (8.20), and is given by 


Po lo fo Ai{[ia„£/«0)]i(r - c,)}dYdY 

[ia o f/o(0)] 5- Ai'{ - [m o t/o(0)] Sc x } 

(8.23) 

Correct matching with the core then requires 


p _ -[«ao^o( 0 )]^oy(0)Ai'{-[zao^o( 0 )]*ci} 

°“ alf™ Ai{[ia o U!>(0)]3{Y -c^dY ’ 

(8.24) 

together with 

h = Pi(0), 

(8.25) 

and 

tVo(V) - YV 0 y(Y)]y^ = 6,(0). 

(8.26) 

It is these conditions that then enable us to determine A(ci) introduced previously, 

^ SS° £ Ai{lia o Ui( 0 )]i(Y - 2 i)dY dY 

[iat/o(0)] 3 AP{-[m o t/o(0)] 3 c a } 

(8.27) 

The procedure is then to determine a solution to (8.11), (8.12) consistent with (8.13), (8.14), 
and (8.27). However, the problem for cj is highly nonlinear, and as such is a difficult 
numerical task. However, if we suppose |ai| — ► 00 , then 

-aifoRsP+dy 

C\ -A A “T Uiot 2 ), 

Jo )[R 1 v+ + R 2 p+]dy 

(8.28) 

which, on account of the nature of p + , u + , i?i, R 2 , R 3 can be shown to be purely real to 
leading order; indeed it is interesting to note that although this sytem does admit complex 
values of C\ for c*i = 0(1), it appears that the above fails to capture the non- neutral nature 
of the inviscid modes for a > qo. 

However, the reason for this is clearly illustrated by considering the nature of the inviscid 
system (5.1) - (5.2) in the limit as a — * a 0 . Suppose that we set 

a = ao + dc, 

(8.29) 

where 

\a\ < |q 0 |, 

(8.30) 

then we expect expansions of the form 


Vo = Vo + avi + a 2 v 2 + 0(a 3 ), 


Po = Po T otpx + a 2 p 2 + 0(a 3 ), 


c 0 = otC\ + a 2 c -2 + +O(o 3 ). 

(8.31) 
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The leading order system is 


Ci{v 0 ,Po } = 0, 
C 2 {v o,po} = 0, 


(8.32) 


(see (8.7), (8.9)), with 

5 o (0) = t> 0 (l) = 0. ( 8 - 33 ) 

At the next order, u a and pi are again determined by means of (8.12), although we shall 
defer discussion of the boundary conditions to be applied to this system, and turn instead 
to consideration of the 0(a 2 ) terms. We find 


. . , UoyVoCl , VoU 0y C2 , UoyVlCi IP\C\ , ,, 2 ] 

Ci{v 2 ,p 2 } = -^j 3 - + —nr~ + ~iJi~ + ~fjr [To ~ M °° ol 


IPO 

^U 2 


= Ra, 


C2 + 


U o 2 

Ci_ 

Vo 


[To - MjJJt) + + 2Mic,p, 


(8.34) 


i rT A 2iao rT a , 2iot$C\ A 2^0 r - - * * 1 

£2 {^ 2 ? P2 } = ~—^o^o ^T-^ 0^1 H ™ ^0 - -«r[— C 1 V 1 - c 2 v 0 ] 

io *0 -*0 -*0 

= Rb* 


(8.35) 


However, the expansions above are not uniformly valid as y — > 0, since U 0 (y) - 0(y) 
in this limit, and hence a breakdown to our approximations must occur, specifically when 

y = O(a). If we set 

y=aY, Y=0( 1), (8-36) 

then wo develops in the following manner (that may be demonstrated to be correct a poste- 
riori) 

t>o = a$o(V r ) + a 2 lna$ 1 (F) + a 2 ^ 2 (F) + --- (8-37) 

It is simple to show that 

$ 0 (F) = -^F, (8.38) 

Cl 

where Ao is some constant; this is consistent with (8.33). We shall defer discussion of $i(F) 
until later, while $ 2 (F) is given by 


$2Y 


W)<&5 


A'iTo(O) 


+ 


I<oYT'( 0) 


[t/'(o)y - ci] [1/5(0 )Y - ci] 2 [c/5(o)F - Ci] 2 T 0 (o)[i/5(o)y - c t ] 


c 2 A o T o (0) 


\UZ{0)Y 2 Ko 


(8.39) 


hmoW-Zi } 2 _ c 1 [c/ 5 (o)y-c 1 ] 2 ’ 


where K i is a constant. 
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Asymptotic analysis of (8.39) reveals that as Y — ► oo 

$2 ~ AoK In Y + AjK + A2 In Y + A3 + • • • , (8.40) 

where the A, are all constants. It is easy to show (see below) that 

A MM ,o A ,s 

0 ua O)To(O) [^(O )] 2 • (8 ' 41j 

However, if $ 2 is to make a proper match with the Y = 0(1) region, we must have 

*1 = A o [^(0)]F. (8.42) 

It turns out that the logarithms are crucial for the generation of non-neutral values of 
Cq. The solution of (8.39) may be written 


$2 = \oK(0)[U'(0)Y-c 1 }£ 


dY 


Mm ~ fii] 


+ 


(8.43) 


= Ao[0o(O)K - c x ]/n 


~ ^o(0)F' 


Cl 


+ 


(8.44) 


where we have only retained terms which are of immediate concern. The crucial significance 
of (8.44) is the presenceof the +*> jump in the value of the logarithm (Mack, 1984, for 
example), as the point Y = ci/£/o(0) is traversed; it is this that generates an imaginary 
component to <J> 2 , given as Y — ► 00 by 


**2~A O 7r[t/'(0)r-Co], 


(here and below a superscript l V denotes an imaginary component). 

Returning to our discussion of the core region, we see we must have that 


(8.45) 


Ui(0) = tfi(l) = 0, (8.46) 

and consequently ci is determined by the expression (8.14), with A = 0. On account of the 
nature of R\ , f? 2 1 R31 v + ,p + we expect c\ to be a real quantity. 

We are now in a position to specify the following boundary conditions 

u 2 (l) = 0, t>«(0) = -AoCitt, (8.47) 

and thus c‘ 2 may be determined, in principle, by means of the condition 

/ [Ra^ + + R bP + ]dy = fi(0)A o c 1 7T, (8.48) 

J 0 

where superscript ‘r’ denotes a real part. 

Unfortunately, we see from (8.41), that A 0 is dependent upon Tq( 0) and Uq (0), and since 
both these quantities are zero in the case of an adiabatic lower wall, then in this case we 
must consider the expressions (8.31) up to 0(a 3 ) in order to determine a value for the 
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imaginary component of the wavespeed as ot * 0; however, since the general technique is 
well established above, we do not carry this out in this paper. 

The key result, therefore, is then that /m{co} is 0(a 2 ) in general (but 0(d 3 ) in the 
insulated wall case), and as a consequence we do not expect to obtain non-neutral values of 
ci, from (8.14) as a x -> oo (this would only be expected in a study of the 0(e 2 ) terms of 
the c expansion, in general, and the 0(e 3 ) terms in the case of an insulated lower wall). The 
above also clearly illustrates how c 0 is only complex for values of a > a 0 , since for a < a 0 
the ‘jump’ in the value of the logarithm is not present. 


9* A Change of Boundary Conditions 


It is curious that the (confined) flow profile under consideration in this paper has such 
different stability characteristics from those of the (external) boundary layer. This leads us 
to question the nature of these fundamental differences. 

One obvious candidate for investigation is the effect of the boundedness of the domain 
(i.e., the impermeability condition imposed on two walls). In an attempt to assess this effect, 
we considered the problem where the impermeability constraint on the upper wall (y = 1), is 
replaced by one of radiation, while retaining the same basic profile. Specifically, we replace 

( 5 -3) by , ni , 

uo(0 = 0 , (9-1) 


(as before), but with 

^ + i/5, o = 0 on y = 1, (9-2) 

oy 

where 

V = ±a[l - Ml{l - Co) 2 }* ■ (9-3) 

Here the sign of v is chosen to ensure that Rl{v} > 0. Indeed this amounts to considering 
the piecewise continuous flow which consists of our basic Couette flow for 0 ^ y 1, and 
uniform flow (U 0 (y) = T 0 {y) = 1) over 1 < y < oo. This may be shown using the arguments 
used by Drazin and Reid (1981) who considered incompressible piecewise linear velocity 
profiles and imposed continuity of pressure at y = 1. 

The system defined by (5.1) - (5.2), (9.1) - (9.2) is solved numerically by the techniques 
used in obtaining the results of Section 6 (i.e., a Runge-Kutta shooting scheme in conjunction 

with a Newton iteration scheme to iterate on Co). 

Results for one example were obtained, namely for M = 2, with adiabatic conditions on 
the temperature at y = 0. Results for Rl{co) are shown in Figure 23. Two modes were found, 
the first (mode I) originates at a = 0, with a wavespeed Co = 1 - 1/M^ = 0.5, and as such 
is typical of so-called “first modes” observed in supersonic boundary layer stability studies 
(e.g., Mack, 1987). Mode II, on the other hand, originates with a wavespeed cq « 1.48, 
and consequently is more typical of the bounded flow modes observed in previous sections 
of this paper. Figure 24 shows the variation of Im {co} with a. Mode I is unstable over the 
entire range of a. shown, while mode II is neutrally stable up to a « 3.5, at which point 
Rl{c 0 } drops below unity, and 7m{co} becomes negative, indicating a stable mode. We 
were unable to find other modes over the range of a considered, particularly modes with 
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m°o) < 0- Thus it appears that the imposition of radiation type boundary conditions 
results in a hybrid situation, with stability characteristics similar to those found in both 
bounded and unbounded flows. 

10. Conclusions 

In this paper, we considered the linear stability of compressible plane Couette flow. Our 
numerical results of the full governing equations are clearly in agreement with the predictions 
of Morawetz (1952, 1954), and also there is a qualitative agreement with the work of Glatzel 
(1989), although the latter author considered a model basic flow. 

The main thrust of this paper has been directed towards modes which have been pre- 
dorninantly inviscid in nature. It appears that these are the most important, since they may 
become unstable, while our studies indicate that viscosity plays a generally stabilizing role, 
throughout. 

The expansions of the assumed form (4.11) are typical of solutions which are predom- 
inantly inviscid, but have viscous corrections (see also Morawetz 1952). In Section 8 we 
studied the important regimes where cq -► 0 (or Co -+ 1) in which the first two terms in 
the series become comparable, a critical layer forms close to one of the walls, and the vis- 
cous layer thickens to 0(Re - ^). These are the “exceptional” cases referred to by Morawetz 
(1952). 

In Section 9, it was shown how, to a large extent, it is the boundary conditions imposed 
on disturbances, that determines the nature of the stability of the flow. Indeed, earlier in the 
paper, in Section 4 it was shown in the study of the G.I.P. condition that it is important to 
impose correct (i.e., radiation type) boundary conditions in the case of unbounded flows to 
avoid the (erroneous) requirement that a G.I.P. is necessary for the existence of supersonic 
neutral modes. 
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Figure 1: Spectrum at M ^ — 2, Re = 2 x 10 5 , at 100 pts. (a) a 
a = 0.20 and (d) blowup of Figure lb. 
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Figure 6: Phase velocity spectrum at = 5, a = 0.1, Re = 5 x 10 6 . (a) 100 pts, (b) 125 
pts, (c) 150 pts. 
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Figure 9: C{\[Rt as a function of a for = 2 (mode II), and the two Reynolds numbers 
j^ e _ 7 x 10 s and Re= = 1.4 x 10 6 . The third curve is Im{ci} obtained from asymptotic 

theory. 



Figure 10: c r as a function of a for M m = 2 (mode II), and the two Reynolds numbers Re 
= 7 x 10 5 and Re= = 1.4 x 10 6 . 
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/m{co} x 10 4 


Figure 12: Im{co } as a function of a for Moq = 2, adiabatic lower wall, mode I. 
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Figure 13. 7m{co} as a function of a for M — 2, adiabatic lower wall, mode II. 


40 


Rt{c 0 } 



41 



/W{co} 



42 





Figure 16: Im{co} as a function of a for Moo = 5, adiabatic lower wall, modes I A , I B , Ic- 
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Figure 17: Im{co} as a function of a for Moo = 5, adiabatic lower wall, mode II. 
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Figure 19: Im{co } as a function of a for M <*> = 5, T( 0) = 1, modes I and II. 


46 














49 



5 1 


NASA 


1. Report No 

NASA CR- 1876 10 


ICASE Report No. 91-64 


4 Title and Subtitle 


Report Documentation Page 


2 Government Accession No 


ON THE LINEAR STABILITY OF COMPRESSIBLE PLANE COUETTE 
FLOW 


3 Recipient s Catalog No. 


5 Report Date 


July 1991 

6 Performing Organization Code 


; 7 Authorise 

Peter W. Duck 
i Gordon Erlebacher 
M. Yousuff Hussaini 

9 Performing Organization Name and Address 

Institute for Computer Applications in Science 
and Engineering 

Mail Stop 132C, NASA Langley Research Center 
Hampton, V A 23665-5223 

12 Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23665-5225 


8. Performing Organization Report No 

91-64 


10 Work Unit No 


- 505-90-52-01 

11. Contract or Grant No 

NAS1-18605 

13. Type of Report and Period Covered 

Contractor Report 

14. Sponsoring Agency Code 


15 Supplementary Notes 

Langley Technical Monitor: 
Michael F. Card 


Submitted to Journal of Fluid 
Mechanics 


Final Repor t 

16 Abstract ~ ~ ~ ~ ~ " — — — 

The linear stability of compressible plane Couette flow is investigated. 

The correct and proper basic velocity and temperature distributions are perturbed 
by a small amplitude normal mode disturbance. The full small amplitude disturbance 
equations are solved numerically at finite Reynolds numbers, and the inviscid limit 
of these equations is then investigated in some detail. It is found that instabil- 
ity can occur, although the stability characteristics of the flow are quite dif- 
ferent from unbounded flows. The effects of viscosity are also calculated, asymp- 
totically, and shown to have a stabilizing role in all the cases investigated. 
Exceptional regimes to the problem occur when the wavespeed of the disturbances ap- 
proaches the velocity of either of the walls, and these regimes are also analyzed in 
some detail. Finally, the effect of imposing radiation-type boundary conditions on 
the upper (moving) wall (in place of impermeability) is investigated, and shown to 
yield results common to both bounded and unbounded flows. 


17 Key Words (Suggested by Authorlsl) 


18 Distribution Statement 


stability, Couette flow 


| 02 - Aerodynamics 


19 Security Ciassif (of this report) 

Unclassified 


r L Unclassified - Unlimited 

20 Security Ciassif (of this page) 21 No of pages , 22 Price 

Unclassified 52 I A q^ 


NASA FORM 1626 OZ y 86 


NASA-Langley, 1991 



